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Abstract We consider the semilinear stochastic heat equation perturbed by additive 
noise. After time-discretization by Euler's method the equation is split into a linear 
stochastic equation and a non-linear random evolution equation. The linear stochas- 
tic equation is discretized in space by a non-adaptive wavelet-Galerkin method. This 
equation is solved first and its solution is substituted into the nonlinear random evo- 
lution equation, which is solved by an adaptive wavelet method. We provide mean 
square estimates for the overall error. 



1 Introduction 

We consider the following semilinear parabolic problem with additive noise, 
du-V -(KVu)dt = f(u)dt + dW, i£^,(e (0,7"), 

u = o, xed9, te(o,T), (l) 

u(-,0) = uo, xgf. 

Here T > 0, 9 C R d , d = 1,2,3, is a convex polygonal domain or a domain with 
smooth boundary d3>, and {W(f)} r >o is an L2(£$)-va\ued Q- Wiener process on a 
filtered probability space (Q,J^,F, {^t}t>o) with respect to the normal filtration 
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{&t}t>0- We use the notation H = L 2 (@), V = {&) with || • || = || • \\ H and (■, •) = 
(•, -)h- Moreover, A: V — > V 1 denotes the linear elliptic operator Au = — V • (kVu) 
for u £ V where k(x) > Kq > is smooth. As usual we consider the bilinear form 
a: V x V — > R defined by a(u, v) = (Au,v) for i/,v£V, and (-, •) denotes the duality 
pairing of V' and V. We denote by e~ tA the analytic semigroup in H generated by 
the realization of -A in H with D(A) = H 2 (3>) f\H^ {9). Finally, f:H—>Hisa 
nonlinear function, which is assumed to be globally Lipschitz continuous, i.e., there 
exists a constant Lf such that 

||/(«)-/(v)|| <L f \\u-v\\, u,v£H. (2) 

In a standard way we introduce the fractional powers A@ , j3 eR, and the corre- 
sponding spaces = D(A$) with norms j|v \\p = |A^v||. We also use the spaces 

L2(Q,IiP) with the mean square norms HvH^^^j = (E[||v||^]) 1 . 
We assume for some j3 > that 

IIa^G^Ihs <~, u eL 2 {n,H p ). (3) 

Here Q is the co variance operator of W and | [ • 1 1 hs denotes the Hilbert-Schmidt norm. 

It is known ||8) that (Q3 has a unique mild solution, which is defined to be the 
solution of the fixed point equation 

u(t)=e- ,A u + [' e- { '- s]A f(it(s))ds+ [' e^'-^ dW(s). (4) 
Jo Jo 

This naturally splits the solution as u = v + w, where w is a stochastic convolution, 

w(t)= [' e- {r - s)A dW(s), (5) 
Jo 

which is the solution of 

dw +Awdt = dW, < f < T; w(0) = 0, (6) 
and v is the solution of the random evolution equation 

v+Av=f(v + w), 0<t<T; v(0)=u . (7) 

Our approach will be to first compute w and then to insert it into © which we 
then solve for v. Finally, u = v + w. For the numerical solution we use Rothe's 
method, where we first discretize with respect to time and then discretize the re- 
sulting elliptic problems with wavelet methods. 

Thus, we fix a time step t > 0, set t„ :~ nx with = T, and consider a backward 
Euler discretization of (Q])- With u" mu(t n ) and increments AW" = W(t„) — W(f„_i) 
this reads 

u" + TAu" = u"~ l + Tf{u")+AW", l<n<N; u° = u . (8) 
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Then we decompose u" = v" + w" to get time-discrete versions of © and Q: 

w" + tAw" = w"- 1 + AW", l<n<N; w° = 0, (9a) 

v" + xAv" = v"- 1 + t/(v" +w n ), l<n<N; v° = u . (9b) 



This allows us to solve the linear problem d9"at first and use the result as an input for 
the nonlinear problem j9bl . Moreover, the stochastic influence in j9bl is smoother 
than in d9al ), which allows us to use fast nonlinear solvers. 

We consider now the spatial discretization of ((9). To this end, let Sj be a mul- 
tiresolution space of order m and let {w'}}^ =0 C Sj be the corresponding Galerkin 
approximation of {w"}^ =0 , i.e., 

w'} + xA]w"j = w n f l +PjAW\ \<n<N; w$ = 0. (10) 

We refer to Subsect. [4.3l for details. We enter this approximation instead of w" into 
i9b\ . The corresponding equation reads 

v" + xAv" = v"- 1 + T/(v" +w n j), 1 < n < N; v° = w . (11) 

For each (0 € £2 and for each n > 1 the nonlinear equation in ( fTTT i is solved by an 
adaptive wavelet algorithm to yield an approximate solution v" with tolerance e„: 

||v?-v"|| <e n , 0<n<N. 

Such an algorithm is described in Subsect. l-Ol Note that we use the same tolerance 
for each ft). The output of the computation will then be the sequence 

u n e =v n e +w}, 0<n<N. (12) 

The total error is u\ - u(t„) = (v" £ - v") + (v" - v") + (wj - w n ) + (u n - u(t„)). The 
contributions are bounded as follows, where the constants C depend on | j uq 1 1 Lj , n ^ , 

||At~2i ||hs, and T, referring to assumption (0. 

In Sect. [2] Theorem |2.5l we study the time-discretization error and show 

m^Ju"-u(t n )\\^ (n , H) <C^>l\ (13) 

In Sect. [3]we analyze the perturbation of the nonlinear term: 

„™ II v" ~ AW{D0) < C™* T H - w"\\L 2[ n,Hy d4) 

The multiresolution approximation of the time-discrete stochastic convolution is 
studied in Sect. \5\ Theorem l5.3l 

max W}-w"\\ L2{n , H) < C 2- J ^- m l (15) 
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Finally, for the adaptive wavelet approximation we have 




Therefore, our main result is the following. 

Theorem 1.1. Assume (0 for some j3 > 0. Let {wj}^ =0 C Sj be computed by a 
multiresolution Galerkin method of order m and {v^}^_ by an adaptive wavelet 
method with tolerances £„. Then the total error in (112b is bounded by 



The literature on numerics for nonlinear stochastic parabolic problems is now 
rather large. We mention, for example, 11311 on pure time-discretization and lfTTl[T6l 
on complete discretization based on the method of lines, where the spatial discretiza- 
tion is first performed by finite elements and the resulting finite-dimensional evolu- 
tion problem is then discretized. Wavelets have been used in IflOl where the spatial 
approximation (without adaptivity) of stochastic convolutions were considered. 

Our present paper is a first attempt towards spatial adaptivity by using Rothe's 
method together with known adaptive wavelet methods for solving the resulting 
nonlinear elliptic problems. 

The spatial Besov regularity of solutions of stochastic PDEs is investigated in 
0[3|. The comparison of the Sobolev and Besov regularity is indicative of whether 
adaptivity is advantageous. For problems with additive noise on domains with 
smooth or convex polygonal boundary, where the solution can be split as u = v + w, 
we expect that the adaptivity is not needed for the stochastic convolution w. We 
therefore apply it only to the random evolution problem ©. Once the domain is not 
convex or the boundary is not regular adaptivity might pay off also for the solution 
of the linear problem (l9al >. 

The recent paper (TJ is a first attempt to a more complete error analysis of Rothe's 
method for both deterministic and stochastic evolution problems. The overlap with 
our present work is not not too large, since we take advantage of special features of 
equations with additive noise. 



2 Pure time discretization 

In the proofs below we will often make use of the following well-known facts about 
the analytic semigroup e~ fA , namely 



max \\ul-u{t„)\\ L2(nM) 



< CT min ^4> +C 2- Jmia{fijn) +maxe, 



n 



where C = C(\\u Q \\ L2{nfilj) 



! ||A^<2^||HS,n 



110 



e 



tA 



\A a &-' A || <Cr a , a>0,t>0, 
-I)A- a \\<Ct a , 0<a<l,r>0, 



(16) 
(17) 
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for some C — C(a), see, for example, El Chapter II, Theorem 6.4]. Also, by a 
simple energy argument we may prove 

f \\A?e- sA v\\ 2 ds<hv\\ 2 , veH,t>0. (18) 
Jo 2 

We quote the following existence, uniqueness and stability result from |9| Lemma 
3.1]. For the mild, and other solution concepts we refer to El Chapters 6 and 7]. 

g-i i 

Lemma 2.1. If ||A^~g 2 || H s < °°for some j3 > 0, mo £ Lo(Q,H), and yji holds, 
then there is a unique mild solution {u(t)} t >o of ([T} with sup ;e [ 0r ]E||M(f)|| 2 < K, 
where K = K(uq, T,Lf). 

Concerning the temporal regularity of the stochastic convolution we have the 
following theorem. 

Theorem 2.2. Let \\A^^Q^- ||hs < °° far some rj £ [0, j]. Then the stochastic convo- 
lution w(t) : = Joe- ( '- a ' A dW((7) is mean square Holder continuous on [0, °°) with 
Holder constant C ~C(rj) and Holder exponent ^ — 77, i.e., 

(E\\w(t)-w(s)\\ 2 y <C\t-s\^, t,s>0. 

Proof. For 77 = 5 the result follows from Lemma |2~T1 Let 77 £ [0, 5) and, without 
loss of generality, let s < t. By independence of the increments of W, 



E\\w{t)-w{s)\\ 2 = E J i 



-(t-a)A 



dW(a)\ 



e -(*-a)A_ e -( s -a)A dw((y )|| 2 =/ 1+ / 2 . 



From Ito's isometry and ([T6l l it follows that 



< cj\t a)- 2 " ||A-"gi da < T -^(f - sf- 2 ^ \\A-"Ql \\i s - 
Finally, let {e^j^Lj be an orthonormal basis of H. Then, by (fTTl i and (TlSTt . 

h = 

= E I" \\(e- {, - s)A -I)A-^-^A^e-^ A Qie k \\ 2 da 
<C(f-5) 1 - 2? '£ r||A2e-(' s - <7 > A A- , 'e2 eit || 2 dc7 



(e -(r-a)A_ e -( J -<TM )e7 ||2 isd(J 



k=\ 



(19) 



< C(t - s) 1 - 2 * £ ||A-"<2^*|| = C(t - s) 1 - 2 * HA^fis I 
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The next result shows that the time regularity of w transfers to the solution of the 
semilinear problem. 



Theorem 2.3. If u g L 2 (Q,H^) and ||A Vg3 || HS < °o/ or some < J3 < 1, then 
there is C = C(T, uq,[3) such that the mild solution u of ^ satisfies 

(E\\u(t) - u(s)\\ 2 y <C\t-s\^, t,se [0,7]. 
Proof Let T > and < s < t < T. Then, by ©, 

u{t)-u{s) = (e-' A -e- sA )u + J' e-('-'^ A f(u(r))dr 

+ /V 
Jq 



*-{l-r)A_-{s- 



r ^)f(u(r))dr + W (t)- W (s). 



In a standard way, for < j3 < 2, we have E||(e-' A -e- iA )«o|| 2 < C\t -i|%||«o|||. 
Using that / is Lipschitz and hence ||/(m)|| < C(l + ||w||), the norm boundedness of 
the semigroup e~' A , and Lemma lXTl we have that 



E 



f e -('- r ) A f(u(r))dr 2 <C|f-^| 2 (l+ sup E|| M (r)|| 2 ) <C\t-s 
Js V re[0,T] 



For < J3 < 1, by Lemma|TT] ^ and (Qj]i, it follows that 



E 



(e 



<sE 



-(t-r)A _ -{s- 



r ^ A )f(u(r))dr 



(e -(-'-M_e-('-'-M) / ( M (r))|| 2 dr 



<Cs(l+ sup E\\u(r)\\ 2 ) \\e- { '- r)A -e- {s - r)A \\ 2 dr 

V re[0,T] ' J 

< Cs f \\ A h-^ A (e-^ A -7)A-! || 2 dr < C\t-sfs 2 -P <C\t-sf. 



Finally, by Theorem |2~21 with 77 = we have E\\w{t) - w(s) || 2 < C\t - s\& , 

which finishes the proof. 

In order to analyze the order of the backward Euler time-stepping ^ we quote the 
following deterministic error estimates, where r{xA) = (1+ tA) -1 . 

Lemma 2.4. The following error estimates hold for t n = nx > 0. 



-nzA n 



r"(TA)]v||<CTT||v||p, 0<j8 < 2, 



-nzA ji 1 



r n (TA)]v\\<CTt- l \\v\l 
£t [/(TA)-e-^]v ^Ct^IIvH 2 ^, < J3 < 2. 



(20) 
(21) 

(22) 



k=\ 
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Proof. Estimates < f20b and ( |2TT i are shown in, for example, 0141 Chapter 7]. Estimate 
d22l ) can be proved in a similar way as (2.17) in |fl6l Lemma 2.8]. 

Theorem 2.5. If u 6 L 2 (Q 1 H 13 ) and || H s < 00 far some < j3 < 1, then 

there is C = C(T,uq,[5) such that for < T < j 2 -, the solutions u of © and u" of 
(O satisfy 

(E||M(f„)- M "|| 2 )2 <Ct^ 2 , t n =nTE[0,T}. 
Proof. We have, with e" := u(t„) — u n , 

e" = [e-'" A - r n {TA)] U Q + £ f * [ e - i! "-^ A - r"- k+l (tA)] dW{s) 

k=f' t i-l 

+ t r ^ {, "- s)A f(u{s))~r"- k+ \TA)f{u k )&s = e l+ e 2 + e,. 

k=\ J 'k-l 

The error e\ is easily bounded, using (|20| |. as 

E|ki|| 2 < Ct^E||m ||^, 0<j3<2. 
The contribution of e 2 is the linear stochastic error. First, we decompose e 2 as 

e 2 = L I'" [C'- k+lA -r"- k+l (TA)]dW{s) 

+ L T [e-('»- s ) A - e -'-*+' A ] dW(j) = e 21 + e 22 . 

k=l Jt k-l 

Let {ff\Y=\ be an ONB of H. By Ito's isometry, the independence of the increments 
of Wand (23, 

E|k2i|| 2 = f T||[/*(TA) -e-^Hls <= £ £ r||[r*(TA) -e-^](2V/|| 2 
*=1 i=U=l 

< cfy iieV/ii 2 ,-! = ct" iia^g^ His, o < p < 2. 
i=i 

The term e 22 can be bounded using a similar argument as in (T% by 

E|| e22 || 2 <CT0||A^<22||HS> 0<j8 <2. 
Next, we can further decompose ej, as 
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e 3 = E r n ~ k+1 (TA) \f(u(t k ))-f(u k )} 6s 

k=\ J 'k-l 

+ t I'" [e-''- t+lA -r"-' ;+1 (TA)]/( M (^))d5 

k=\ J 'k-l 

+ t /'* e-<>-^ A lf(u(s))~f(u(t k ))]ds 

k=l J 't-l 

+ L /'* [e" (r "^ )A -e- r "-* +lA ]/(«(^))d^ fi31 + e 3 2 + e 3 3 + e 3 4 . 
By the stability of r"(rA) and the Lipschitz condition on /, we have 
E||e 3 i|| 2 < 2L 2 f T 2 E\\e n \\ 2 + 2L 2 f T: 2 n"Y i E\\e k \\ 2 < 2l} f t 2 E\\e"\\ 2 + CtY i E\\e k \\ 2 . 

k=l k=l 

By (EB and Lemma [27TI with Tt~} k+1 = {n - k+ 1) -1 = 

E|| e3 2|| 2 < CE( £ t«-^ +i ||/("(f,))ll) 2 < Ct 2 £ i £ E||/(«( f ,))|| 2 

n 

< Ct 2 +E||«(f it )|| 2 ) < Crf„ < Ct. 
jt=i 

Furthermore, by Theorem l2.3l 



E|k33|| 2 <f„E / nfHs))-f(u(tk))\\ 2 ds<Ct 2 ^ <CtP,0<I3<1. 
To estimate we have, using again that t n _k+\ =t n — f>t l an d Lemma [2~Tl 



2 



E|k 3 4|| 2 = E(£ f !|[A^e-('«- s ) A (/-e-( s -'*-i) A )]A-4 / ( M ( J ))||d,) 

<a„f r (^-^-^^eh/^))!! 2 ^^^, o</3<i. 
k=.\ J h-\ 

Putting the pieces together, we have 

E||el 2 <C^+24T 2 E|| e "|| 2 + CT£ E l! ei H 2 > < /3 < 1 . 

>t=i 

Finally, if T < 277 > then by the discrete Gronwall lemma, 

E||e„|| 2 < CT P e Ct " < Ct P , < j3 < 1, 
and the theorem is established. 
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3 Error analysis for the nonlinear random problem 

In this section we bound the term E[||v" — v"|| 2 ] in (fl4] >. We use the global Lipschitz 
condition (0. 

Lemma 3.1. Assume that %Lf < A. Then, with C = 2LfTe 2L f T , 



max 

\<n<N 




Proof. Let e n :=v"-v". Then, we have by (gb]) and dTTb 

e" + TAe" = T (/(v»+wJf)-/(v" + w")) +e n ~ 1 . 
Since e° = 0, we get by induction 

e" = X £ (/ + tA)-(" +1 -^ + wj) - /(W + w } )) . 
In view of the global Lipschitz condition (O, this results in the estimate 

||*"|| < LfX £ || (7 + TA)-(" +1 -» II || v> + W } - v> - w>|| 
;=i 

<^t(ik-^'ii+iki), 

;=i 

since ||(/+tA) _1 || < 1. Thus, we obtain 

ii e »n < (i -7/t)- v( £ ik - +"f iki) . 

7=1 7=1 

With L/T < i we complete the proof by the standard discrete Gronwall lemma. 



4 Wavelets 

4.1 Wavelet bases 

In this section, we collect the notation and the main properties of wavelets that will 
be needed in the sequel. We refer to J4] [7] [T3] for more details on wavelet methods 
for PDEs. For the space discretization, let 
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be a biorthogonal basis of H, i.e., in particular (ft t Y^ H = dg^- Here, X typically 
is an index vector X = (j,k) containing both the information on the level j = \X\ 
and the location in space k (e.g., the center of the support of ft). Note that f also 
contains the scaling functions on the coarsest level that are not wavelets. We will 
refer to \X \ — as the level of the scaling functions. 

In addition, we assume that ft G V, which is an assumption on the regularity 
(and boundary conditions) of the primal wavelets. To be precise, we pose the fol- 
lowing assumptions on the wavelet bases: 

1. Regularity: ft G H'(@), X G J w for all < t < s w ; 

2. Vanishing moments: ((•)', Wx)o-,S> = 0, < r < m*y, \X\ > 0. 

3. Locality: diam(supp ft) ~ 2 - ' L 

We assume the same properties for the dual wavelet basis with s»p and mvy replaced 
by sip and m<p. Note that the dual wavelet ft does not need to be in V, typically one 
expects ft e F'. 

We will consider (often finite-dimensional) subspaces generated by (adaptively 
generated finite) sets of indices A C ^ and 

f A := {ft : X G A}, S A := c^pan^), 

where the closure is of course not needed if A is a finite set. If A = Aj := {X £ 
J?^ ' '■ \X \ < /— 1}, then Sj := S A , contains all wavelets up to level J — 1 so that Sj 
coincides with the multiresolution space generated by all scaling functions on level 
J, i.e., 

S^span^, <Pj = {(pj^ : k G Jj), (23) 
where J>j is an appropriate index set. 



4.2 Adaptive wavelet methods for nonlinear variational problems 

In this section, we quote from J6| the main facts on adaptive wavelet methods for 
solving stationary nonlinear variational problems. Note, that all what is said in this 
section is taken from j6] . However, we abandon further reference for easier reading. 

Let F : V — >• V 1 be a nonlinear map. We consider the problem of finding u G V 
such that 

(v,R(u)):=(v,F(u)-g)=0, veV, (24) 

where g G V 1 is given. As an example, let F be given as (v,F(u)) := a(v,u) + 
(v,f(u)) which covers ( fTTT ). The main idea is to consider an equivalent formula- 
tion of d24l i in terms of the wavelet coefficients u of the unknown solution u = u r f . 
Setting 

R(v):=((ft, J R(v))) AG/V , v = v r f, 
the equivalent formulation amounts to finding ue^t/ 5 ') sucn mat 
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R(v)=0. (25) 

The next ingredient is a basic iteration in the (infinite-dimensional) space lii^^) 
and replacing the infinite operator applications in an adaptive way by finite approx- 
imations in order to obtain a computable version. Starting by some finite the 
iteration reads 

u («+i) =u W_ 4u W 5 Au (n) := bWr( u W) (26) 

where the operator fiW is to be chosen and determines the nonlinear solution 
method (such as Richardson or Newton). The sequence AvS"' = BWR(uW) (pos- 
sibly infinite even for finite input u™) is then replaced by some finite sequence 

:=RES[j] n ,BW,R,uW] such that 

||V")-w< n) ||<77,, 

Replacing AvS"^ by wl in d26l ) and choosing the sequence of tolerances (J7,,);igN 
appropriately results in a convergent algorithm such that any tolerance e is reached 
after finitely many steps. We set u(e) := SOLVE[£,R,BM,u( '] such that we get 
|u-u(e)|| < e. 

In terms of optimality, there are several issues to be considered: 

• How many iterations n(e) are required in order to achieve £-accuracy? 

• How many "active" coefficients are needed to represent the numerical approxi- 
mation and how does that compare with a "best" approximation? 

• How many operations (arithmetic, storage) and how much storage is required? 

In order to quantify that, one considers so-called approximation classes 

/:={v64(/):Mv)<r'} 

of all those sequences whose error of best N-term approximation 

C7jv(v) := min{||v — w||f 2 : # suppw < N} 

decays at a certain rate (suppv := {A £ : ^ 0}, v = (vx)xe J^- ^ or 
case V = H' it is known that the solution being in some Besov space u G Z? r T +</s (L T ), 
T = (i+ i) — ^ implies that u G si s . 

Let us first consider the case where F = A is a linear elliptic partial differential 
operator, i.e., An = g £ V', where A : V — > V' , g £ V 1 is given and u £ V is to be 
determined. For the discretization we use a wavelet basis f in H where rescaled 
versions admit Riesz bases in V and V', respectively. Then, the operator equation 
can equivalently be written as 
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where A := D" 1 a( , P,»F)D _1 , g := D^ 1 (g,'F) and u := D(ux) Xe *v, with u x being 
the wavelet coefficients of the unknown function u S V, \\u\\y ~ ll u ll^ 2 ( a / s ')' Wavelet 
preconditioning results in the fact that K"2(A) < °°, 0. 

The (biinfinite) matrix A is said to be s* -compressible, A € if for any < 
s < s* and every j 6 N there exists a matrix A ; - with the following properties: For 
some summable sequence (a,) 7£ N, the matrix Aj is obtained by replacing all but 
the order of up.' entries per row and column in A by zero and satisfies 

||A-A ; -]|<Ca;2-^ jeN. 

Wavelet representations of differential (and certain integral) operators fall into this 
category. Typically, s* depends on the regularity and the order of vanishing moments 
of the wavelets. Then, one can construct a linear counterpart RESi; n of RES such 
that Wr) := RESr m [?7,A,g,v] for finite input v satisfies 

K-(Av-g)|| f2 < 77, (27a) 
KlU/^CIvll^ + llulU), (27b) 
#suppw,<77-^(||v||^ + ||u||^), (27c) 

where the constants in ( 127bl i. ( I27cb depend only on s. Here, we have used the quasi- 
norm 

||v||. c/s := s\xpN s a N (\). 

NeN 

This is the main ingredient for proving optimality of the scheme in the following 
sense. 

Theorem 4.1 (|5, 6|). If A E and if the exact solution u of Au = g satifi.es u € 



gf s , s < s*, then u(e) = SOLVE [in [e] satisfies 

||u-u(e)|| <£, (28a) 

#suppu(e) < g- 1/s , (28b) 

computational complexity ~ #suppu(e). (28c) 



It turns out that most of what is said before also holds for the nonlinear case 
except that the analysis of the approximate evaluation of nonlinear expressions R(v) 
poses a constraint on the structure of the active coefficients, namely that it has tree 
structure. In order to define this, one uses the notation ji -< X, X,jX £ ^^to express 
that /J is a descendent of A. In the univariate case with i//^ = y/j ^ = yl 2 y(V .-k) 
this is obvious since (j + 1, 2k) and (j + 1,2A:+ 1) are children of X. In higher 
dimensions and even on more complex domains this can also be defined - with 
some more technical effort, however. 

Then, a set ^ C is called a tree if A £ 2? implies [I e 2? for all jx <E 
with X < ji. Given this, the error of the best N -term tree approximation is given as 

°iv ee ( v ) : = min{||v- w 1 1 f 2 : 5* :=#suppw is a tree and#^ < N} 
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and define the tree approximation space as 

<e:={ve4(/):C(v)<r 5 } 

which is a quasi-normed space under the quasi-norm 

||v|| <ee :=sup^c7r(v). 

The extension of the s* -compressibility c & s * is the s*-sparsity of the scheme RES 
which is defined by the following property: If the exact solution u of (|23T > is in s^fee 
for somes < s*, then Wjj :~ RES[tj,B,R, v] for finite v satisfies 

K IL< CC < c(||v||, <cc + ||u|| <cc + i), 

#suppw,<C77- 1/s (l|v||^ + ||u||^ +1), 
comp. complexity ~C(rj-^ s (\\v\\ l i' + ||u|| 1 jj + 1) +#^(suppv)), 

tree tree 

where C depends only on s when s — > s* and 5^(suppv) denotes the smallest tree 
containing suppv. Now, we are ready to collect the main result. 

Theorem 4.2 (f6, Theorem 6.1]). //RES is s* -sparse, s* > and ifu€ &?4eef or 
some s < s*, then the approximations u(e) satisfy |u — u(e) || < £ with 

#suppu(e) < Ce-V'IMI^, l|u(e)|U Cc < C||u|| <ec , 
where C depends only on s when s — > s*. The number of operations is bounded by 



4.3 Multiresolution Galerkin methods 

Let Sj = Sj\j be a multiresolution space (f23t . The multiresolution Galerkin approx- 
imation of the equation Am = / in V' is to find uj &Sj such that 

a(uj,vj) = (f,vj) VveSj. (29) 

Define the orthogonal projector Pj\ H —> Sj by 

(Pjv, wj) = (v, wj), v£H,wj£Sj. (30) 

Note that Pj can be extended to V' by ( f30b since Sy c V. Next, we define the operator 
Aj: Sj^Sjby 

a(AjVj,wj)=a(vj,wj), uj,vjESj. 
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Then ( f29T > reads Ajuj — Pjf in Sj. Alternatively we may write uj = Rju, where 
Rj : V — > Sj is the Ritz projector, defined by 

a(Rjv,wj) = a(v,wj), v G V, wj e Sj. 

The multiresolution space is of order m if 

inf ||v-w / ||<2-'" y ||v|| m; ^, veH m (9)nV. 

WjESj 

Standard arguments then show, using elliptic regularity thanks to our assumptions 
on 3$, that \\uj — u\\ < 2~ mJ \\u\\ m -^, or in other words 

||v-^v||<2-' tt/ ||v|| m; ^ vGH m (9)nV. (31) 



5 Error analysis for the stochastic convolution 



The next lemma is of independent interest and we state it in a general form. 

Lemma 5.1. Let —A and —B generate strongly continuous semigroups e~ tA and 
e~' B on a Banach space X and let r(s) = (1 + s) . Then, for all x,y £ X, N € N, 
T>0, 



n=\ 



t£ \\r n {xB)y-r n {xA)x\\ p < / \\e-' B y - e-' A x\\ p dt , 1 < p << 



\r n {%B)y-r n {%A)x\\<su V \\z- ,a y-z- tA x\\. 
t>0 



(32) 
(33) 



Proof. By the Hille-Phillips functional calculus, we have 



r n {TB)y-r n {xA)x= / (e^y ~ e -" A x)f n (t)dt, (34) 
Jo 

where /„ denotes the «th convolution power of /(f) = e - '. Since ||/,,||l 1 (r + ) = 1 
inequality d33l ) follows immediately by Holder's inequality. To see d32b we note that 
/„ is a probability density and hence by Jensen's inequality and ( T341 ). 



T£||r"(TB)y-r"(TA)*r = T£ / {<T nB y -<T nA x)f n {t)dt 



n=\ 



n=\ 



<t£ / \\e-" B y-e-" A x\\Pf n (t)dt 



||e-' B y-e-'^||"df sup £/„(*). 

r>0„=i 



Finally, by monotone convergence, the Laplace transform of £™ =1 /« i s gi yen by 
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{ifn)\v=im=i( T ^y=l a>o. 

n=l n=\ n=l 

Thus, Y,n=i fn = 1 an d the proof is complete. 

Next we derive an error estimate for the multiresolution approximation of the 
semigroup e~ tA and its Euler approximation r"(xA). 

Lemma 5.2. Let Sj be a multiresolution space of order m and let A, Aj, and Pj be 
as in Subsection \4.3\ Then, for T > 0, N > 1, T, we have 

(jf ||e-' A 'P / v-e-' A v|| 2 d f ) 2 <C2-^||v|| /3 _ 1 , 0<j3<m, (35) 

and 

(t£ V{xAj)Pjv- r"(TA)v|| 2 ) 2 < C2-^\\v\\p_ h 0<j3<m. (36) 

Proof. Estimate d35l ) is known in the finite element context, see for example |[T4l 
Theorem 2.5], and may be proved in a completely analogous fashion for using the 
approximation property ( BIT ) of the Ritz projection R j, the parabolic smoothing dTlt . 
and interpolation. Finally, (l36l l follows from d35l l by using Lemma IBTTI with x = v, 
y = Pjv, and B =Aj. (Note that C is independent of T.) 

Now we are ready to consider the multiresolution approximation of w" in d9al >. 

Theorem 5.3. Let Sj be a multiresolution space of order m and w and w'j the solu- 
tions of d9at and ( fTOt . ||AT~g2 ||hs < °°for some < j3 < m, f/zen 

(E[||w}-w"|| 2 ])i<C2-^||A fi T i e2|| HS . 
Proo/ Let ^ = kx, k = 0, ...,n. By JTOl . ( |9al ), and induction, 

- ^" = £ T [f* - *" 1 (TAj)Py - r""* +1 (TA)] dW{s), 

k=\ Jt k-\ 

whence, by Ito's isometry, we get 

nW}-^\\ 2 ] = i T \\[r"- k+l (rAj)Pj-r"- k+ \TA)]Qi\\ 2 m ds 

k=\ Jt k-\ 

= Y.x\\[r k (xA J )P J -/{xA)}Q^ m . 

k=\ 

Let {ei}Y = i be an orthonormal basis of H. Then, using Lemma l5T2l we obtain 
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E[\\w"j-w' 



\ 2 } = HT\\{r k (TAj)Pj-r k (TA)}Q?e l \ 



<C^2-^W ei \\l 



1 — 



C2- 2J P\\A^Q±\\l 



i=\ 
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